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Abstract: MicroPattern Gaseous Detectors (MPGD) made of resistive strips have raised as a 
promising technology for the protection against spark processes having place in the gaseous cham- 
ber. The reproduction of the signals and its propagation through the resistive foil is mandatory to 
better understand its behavior and optimize the key parameters which might depend on the applica- 
tion requirements. In this work it will be presented a resistive-strip model and the charge diffusion 
through the resistive strip for different model parameters, such as the strip linear resistivity and 
capacitance, together with the advantages and/or disadvantages of this type of technology. 
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1. Introduction 

High amplification gains are required in MicroPattern Gaseous Detectors (MPGD) in order to 
achieve optimum signal to noise ratios, increasing the performance of detectors in terms of energy 
resolution and efficiency. The gain applied allows to observe signals from gas ionizing interactions 
which produce few primary electrons at the detector conversion region. The populated electron 
avalanches achieved with these few primary electrons entail the risk to produce a discharge at the 
cathode of the detector when the Raether critical electron density is reached, around ~ 10^ — 10^ 
electrons per avalanche [U]. Once a paiticular detector gain has been fixed, different type of in- 
teractions inside the chamber (or produced in the chamber structure) can produce a spark with a 
certain probability which depends on the number of primary electrons generated [ED- 

Nowadays, gaseous detector technology is limited by the presence of sparks in hazard and/or 
high particles flux environments, due to the increased frequency of appearance. Detector discharges 
might affect the detector response in different ways; reducing its operating lifetime due to intense 
currents produced in short periods of time, heating and melting the materials at the affected regions, 
damaging the read-out electronics which have to support huge current loads in a brief period of 
time, and moreover increasing the detector dead-time given that spark phenomena entail the dis- 
charge of the cathode and therefore the amplification field is lost during a relatively long period 
of time, which is required by the high voltage power supply to restore the charges and recover the 
nominal voltage. 

It was first observed in RPC-type^ detectors that the introduction of a high impedance resistive 
coating at the anode limits the detector cun^ent during a period of at least some /xs, constraining the 
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spark process to the streamer phase and reducing the total amount of charge released (ED- Further- 
more, the limited discharge current affects only locally avoiding the field to be lost at the remaining 
detector regions, and thus reducing the dead-time of the detector. 

Micromegas detectors were introduced already in 1995 [M] as a good candidate for high par- 
ticle flux environments, and spark studies with standard micromegas technology were also carried 
out O- Recently, additional efforts are pushing the development of resistive strip micromegas de- 
tectors in order to increase its robustness and protection in hazard or high flux environments by 
limiting spark discharges in the same way as it was done for RPCs. In particular, the MAMMA^ 
collaboration is developing large ai^ea micromegas detectors, which introduce the resistive coating 
technique 0], for the upcoming upgrade of the HL-LHC^, where the luminosity will be increased 
by a factor 10 reaching up to L = lO^^cm^^s^'. 

The MAMMA collaboration has investigated new detector prototypes, with different resistive 
coating topologies. These studies have increased the robustness and stability of micromegas detec- 
tors in the presence of an intense and highly ionizing environment, limiting the negative effect of 
sparks could have on them 

The present work pretends to give a theoretical description of the signal propagation through 
a resistive strip as it was previously done in RPC-type detectors IS] . A resistive strip model will 
be introduced based on the resistive-strip detectors developed in the framework of the MAMMA 
collaboration. The intention of this work is to provide a model of the electric behavior of the 
detector to a current signal in order to better understand the different responses of the detector 
which could be produced in different environments, and moreover the response dependency on the 
key parameters of the detector, which will be defined in the model. It is expected, that the future 
agreement between model and experimental measurements could allow to exploit the still, if any, 
unexplored capabilities of this new type of technology. 

2. The resistive strip model. 

A simplified approach of a single resistive strip is taken into to account in order to calculate the 
induced signal at the micromegas usual read-out and the charge diffusion through the resistive strip. 
The model here described is based on the equivalent electric circuit introduced at [0]. This equiva- 
lent circuit describes in general a standard micromegas detector covered by a thin insulating layer 
where resistive strips are placed facing the standard metallic strips, creating an inter-capacitance 
and therefore an electric transmission line. 

In order to calculate the charge diffusion, the resistive strip is subdivided in a set of differential 
elements which final element is grounded through a boundary resistor R/,. The resistive strip has a 
linear resistivity Rx and is electrically connected to the standard read-out, defined by Vc, through 
the mentioned insulating layer which introduces the strip linear capacitance Cx ■ The standard read- 
out is connected in this model to a low impedance value Rstnp, and the residual capacitance of the 
PCB board Cp^b (see Fig. [I]). 
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Figure 1. Simplified RC-model of a resistive strip. 



In this simplified model, applying Kirchoff continuity at each node we can describe the resis- 
tive strip potential V{x,t) in terms of the current induced at the inter-capacitance Cx 

V{x + 5x,t)-V{x,t) V{x,t)-V{x-5x,t) fdV{x,t) dVc{t)\ 

-Cxoxl — — (2.1) 



Rxdx Rx8x \ dt dt 

in the case where 5x — )• the above equation is reduced to the following expression 



d^VM diVix,t)-V,it)) 

-^-^ = CxRx (2.2) 

a partial differential equation (PDE= which describes the spatial and temporal evolution of the 
potential in the resistive strip, and thus the charge distribution in the strip which is directly related 
with the potential V(x,f) by the relation 



Qx{x,t)=Cx{V{x,t)-V,{t)) (2.3) 

where Qx is the linear charge density in the resistive strip. 

The evolution of the potential and/or charge in the resistive strip is bounded by the external 
circuit elements attached to the strip electrode, Rstrtp and Cpcb, and the fixed resistive value, Rf,, 
attached to the resistive strip. The first elements (Rstrtp and Cpct) allow to determine the strip 
voltage Vc as a function of the resistive strip potential distribution, by integrating the total equivalent 
current induced at the capacitive surface along the resistive strip length xl, 

'c>. = Cx r ^XVM - Vc{t))dx = CpJ-^ + ^ (2.4) 

JO at at Kstrip 

relation that can be rewritten in the most convenient expression 

= r '-^"^ - 1 . (2.5) 

dt Cpcb+XtCxJo at {XLCx+Cpch)Rstrip 

The resistive value Rh provides the boundary condition on the voltage V„{t) at the beginning 
of the strip, and its value is related with the voltage gradient at the resistive strip end. 
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Rh dV{x,t) 
dx 



(2.6) 



The relation ([2.2| ) describes the potential propagation at the resistive strip. However, in order to 
observe such propagation we need a perturbation which is introduced as a time dependent current 
source spatially distributed and described by the change on linear charge density p{x,t) along the 
resistive strip. This current is introduced inside relation ( |Z.2D by taking into account the expected 
voltage increase at the resistive nodes, 



d^V{x,t) 
dx^ 



^^^^ diV{x,t)-V,{t)) dpix,t) 



dt 



dt 



(2.7) 



the introduced charge density pretends to emulate the effect of charge induced by a physical 
avalanche, which would be the connection between our detector read-out and the physical signal. 

The solution of the coupled relations ( |2.5P and ( |2.7| ) for a given input charge density p{x,t) 
will provide a full description of the signal propagation within the frame of this model. 

3. Discretization and resistive strip potential solution. 

If the resistive strip is subdivided in a total of A'^ differential elements of length 8x, the discrete 
model contains a total of A'^ + 1 resistive potential nodes. The discrete relations of the A'^ — 1 internal 
nodes is obtained directly from the relation (^77|), describing the potential evolution in terms of the 
neighborhood potentials and the equivalent current at the given node, 
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where Ti = R^Cx- The potential at the last node, which defines the non-grounded strip end is given 
by the relation 
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and the first strip element which is constrained with the resistor Rf, leads to the following boundary 
relation 
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This set of A'^+ 1 equations can be written in the matricial form 
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which can be reduced to the following expression that represents a set of A/^ + 1 coupled differential 
equations 



dy 1 ^ I dp dVc, 
dt ix8x^ Cx dt dt 



(3.5) 



Furthermore, the potential Vc at the strip read-out, which is described by the relation ([Z.5P, 
depends on the resistive potential V{x,t). In order to introduce this dependency inside relation (p.SP, 
the relation describing dVc/dt must be discretized as follows. 



dVc 
dt 



■trip 



(3.6) 



where = Cx/iCpcb +xlCx). Introducing the last expression inside relation ( ^3] ) we obtain 

1^1 I dp ^Vc ^ 



dt '^x^x^ ~^ '^k^x^ " Cx dt 



CxRsti 



(3.7) 



ip 



where the matrix is defined as 



—E,5x \ — t,5x 



(3.8) 



Finally, the differential dependency on Vc has been substituted by a linear dependency on Vc, 
which can be more easily introduced in the calculation by imposing some initial conditions. If we 
work out the value of dy/dt we obtain the following final expression for the potential change along 
the resistive strip 



dy _ 1 

dt Xx 5x^ 
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dp ^V, 



dt CxR, 



(3.9) 



strip 



expression which represents a set of A'^ + 1 coupled first order differential equations. In order to 
solve the potential it is required to uncouple this set of equations and apply the standard method 
to solve a set of independent first order differential equations. To achieve that the matrix ^ which 
diagonalizes Jif^^ ^ is introduced in the calculation 



-1 



where A is the resulting diagonal matrix. By multiplying to the left side the relation 
and introducing the identity SiS^^^ the term which contains v becomes diagonal 



(3.10) 
by 
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obtaining a set of differential equations which are uncoupled and where the transformed potential 
u = can be solved for each virtual node. Finally, the set of A'^+ 1 uncoupled differential 

equations can be better described by the following expression 

^ = ^Au + ^ ^vo - - TT^^b (3.12) 

where = S^-^J^'K 

The expression ( [3. 121 ) can be numerically solved for a given time step At by any of the popular 
methods (Euler or Runge-Kutta). The time discretization is now trivial in this relation and the value 
of Vc can be updated at each time iteration by relation ( [3.6| ) by calculating first the transformed 
potentials u and obtaining the real resistive strip potential v = at each time iteration, imposing 
the initial conditions to satisfy Vc{t = 0) = and V{xJ = 0) = 0. 



4. Charge diffusion through the resistive strip. 

The model described, together with the method to calculate the propagation, has been used to 
calculate the charge diffusion along the resistive strip. In order to study such propagation a linear- 
current density must be introduced. The signal is chosen to be Gaussian shaped for simplicity. 
Therefore, at each node j, the linear current density can be described by the following expression 

dt ~lna,a, ^\ lo} ) la} ) ^ ' 

where Xg describes the central linear charge density of the virtual event, x^ the mean event position, 
t,„ the time at which the maximum current is reached and, Ox and Ot the spatial and temporal chai^ge 
diffusion previously produced in the gaseous medium. 

A relatively fast signal is generated at t,„ = 200 ns with Ot = 50 ns at the central position, 
Xm = 5 cm with Ox = 0.5 cm. Figure shows the result of charge diffusion due to such current 
signal through the resistive strip, 10 cm long, at different periods of time. 
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Figure 2. Potential distribution along the resistive strip at different periods of time. On the left, 
first periods of time, on the right plot, snapshots after longer time period has passed. 
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This first calculation allows to observe the charge spread along the resistive strip as a function 
of time, and also the asymmetry effect at the boundary resistor which can be observed at the last 
time steps, which is produced by the electrons drifting towards the ground potential. 

The same calculation was performed for different parameter combinations, including linear 
resistivity and linear capacitance Cp^ . Figure shows the effect of those two parameters on the 
charge diffusion, showing a faster charge spread as the resistivity and/or the capacitance Q 
lowered. 
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Figure 3. On the left, charge distribution along the resistive strip for different resistivity values at 
the same period of time. On the right, the distribution for different values of the linear capacitance. 



4.1 Homogeneous current effect on the resistive strip. 

The longer times required for charge diffusion through the resistive strip, and its obvious depen- 
dency with the charge position thi^ough the time constant T;^ [s-cm^'^], will affect the potential ho- 
mogeneity at the resistive anode in high rate and thus intense radiation environments. 

The resistive-strip model previously described has been used to show the effect of an homoge- 
neous illumination on the resistive potential distribution V{x,t). This is achieved by imposing the 
input signal previously given by relation ^?T] to be uniform and constant in time. Figure ^ shows 
the potential V{x,t) transition along the resistive-strip and the effect of the same input cuiTent as a 
function of Rx and Cx on the final state, emphasizing the higher effect with the resistivity. 




Position [mm] Position |mml 

Figure 4. On the left, transition effect on potential distribution due to the homogeneous illumination 
switched on at the starting reference time. On the right, final-state potential distribution for different 
model parameters. 



-7- 



4.2 Resistive charge diffusion dependence on event position. 

The charge diffusion through the resistive-strip becomes necessarily an effective current at the 
boundary resistor, coming from charges in the resistive material that drift towards the ground. This 
current can be easily obtained from the propagation calculated in the model by representing the 
value of the potential Vo{t) at the boundary resistor as a function of time. The pulse-shape of the 
signals propagating through the resistive-strip will depend on the event position, an effect that will 
not be negligible thanks to the time propagation delay induced by the high strip resistivity. Figure ^ 
shows the current-shape at the boundary resistor for similar characteristics events having place at 
different strip coordinates, being the lower position values closer to the strip boundary resistor. As 
it is observed in these plots, the closer the event is to the ground the faster is the expected signal at 
the boundary resistor. 
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Figure 5. Current at the resistive-strip boundary resistor as a function of time, for events produced 
at different positions. On the left plot, linear scale and on the right plot, logarithmic scale. 



Once the event signal is fixed, the different pulse-shapes can be associated uniquely to a given 
position. This property of resistive strips has already been used for position detection in silicon 
based detectors by using the charge division method [Q], introducing an output signal at each end. 
In our particular case, with one unique end, the position sensitivity could be reached by analyzing 
the pulse time properties (see Fig. ^. 
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Figure 6. Pulse time properties, for a resistivity Rx = lOOkO/mm, as a function of the position. 
The relative pulse amplitude is also represented in an auxiliary linear scale. 
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The signal properties can be therefore used to locahze the position of the event in the resistive 
strip. It is remarkable how the pulse amplitude is weakened as the event is farther away from the 
strip-end where the resistor ground is located. This fact will favor the localization of events which 
are near the ground. But at least, events producing huge amount of charge, as developing streamers, 
will produce signals which could still be observed. 

The temporal pulse shape properties depend on the resistive parameters such as resistivity 7?;^ > 
capacitance Cx and boundary resistor values R},. Figure shows the time delay dependency with 
these parameters as a function of the event position. It is observed how the time delay changes 
are increased with higher resistivities and capacitances, and the effect of the boundary resistor 
value is not affecting considerably. A higher change on time delay with position should provide a 
higher equivalent spatial definition, however a compromise between signal delay and signal inten- 
sity should also be considered. 
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Figure 7. Time signal delay as a function of the event position comparing the effect of different 
model parameter values. 



5. Conclusions and future work 

A method for the description of the slow diffusion signal on resistive strip has been formulated. 
Some results concerning relevant matters of the detector have been presented showing the still pre- 
dicting potential of a simplified model. The fuU mathematical description provided for this model 
could allow future extensions to be implemented, as it could be the case considering resistive-strip 
inhomogeneities, and other complexities that cannot be addressed with commercial or dedicated 
simulation software. Moreover, it could provide the key to connect streamer physics and other re- 
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lated physical phenomena with this kind of detectors, allowing to study the response of the detector 
in a variety of physics scenarios and applications. 

New prototypes, based on the MAMMA project ones, have been recently produced by R. 
Oliveira's team at CERN Micromegas workshop. Three different prototypes have been designed 
which combine different resistive strip and metallic strip widths in order to provide a single detector 
with different resistive and capacitive properties. A special contact has been inserted in these 
prototypes in order to have access to the signal coming out from the resistive strips and at the same 
time allowing to change the value of the boundary resistor (see Fig. 




Figure 8. On the left, from top to bottom, detector Gerber design, a picture of resistive strip 
connectors on detector PCB board, and a microscope picture of the resistive strips with different 
widths on top of standard strips. On the right, schematic showing the transversal section of the 
detector. 

A dedicated electronic read-out is being designed for these prototypes in order to preamplify 
and integrate the resistive signals with the adequate timing. The results obtained could allow to 
verify the resistive- strip model described and/or provide clues to apply any modifications required. 
The measurement of resistive propagation signals in different conditions will be studied for the first 
time in this type of detectors. Furthermore, the possibility to read-out the position by using this 
technique could provide a way to reduce the number of read-out channels maintaining good spatial 
detection resolution. 
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